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We investigate theoretically the phenomenon of exchange narrowing in the absorption spectrum 
of a chain of monomers, which are coupled via resonant dipole-dipole interaction. The individual 
(uncoupled) monomers exhibit a broad absorption line shape due to the coupling to an environ- 
ment consisting of a continuum of vibrational modes. Upon increasing the interaction between the 
monomers, the absorption spectrum of the chain narrows. For a non-Markovian environment with 
a Lorentzian spectral density, we find a narrowing of the peak width (full width at half maximum 
(FWHM)) by a factor 1/N, where N is the number of monomers. This is much stronger than 
the usual 1/y/N narrowing. Furthermore it turns out that for a Markovian environment no ex- 
change narrowing at all occurs. The relation of different measures of the width (FWHM, standard 
deviation) is discussed. 



I. INTRODUCTION 

Due to interactions with an environment the absorp- 
tion lines of electronic transitions can be considerably 
broadened. This broadening might be due to the fact 
that the individual (two-level) systems see different envi- 
ronments (static disorder) or that the excitation couples 
to time-dependent fluctuations or vibrational modes of 
the environment (dynamic disorder). In this work we 
refer to the individual two-level systems as monomers. 

When there is exchange of excitation between the 
monomers (e.g. due to transition dipole-dipole interac- 
tion) , excitonic states are formed. In the case of vanishing 
coupling to the environment, this leads to the formation 
of a band of eigenstates having a width proportional to 
the interaction strength between the monomers. Due to 
selection rules not all of the states absorb. For a parallel 
arrangement of the transition dipolcs of the monomers, 
the dominant absorption will be into states at the band 
edge. 

One might expect that for the coupled monomers the 
interaction with the environment will lead to a similar 
broadening as in the case of uncoupled monomers. How- 
ever, for strong coupling between the monomers the ag- 
gregate absorption line is considerably narrowed com- 
pared to that of the uncoupled monomers. This ex- 
change narrowing effect occurs for example in param- 
agnetic resonance in the excitonic line shapes of dis- 
ordered solids 0, or for the J-band of molecular aggre- 
gates [3-6]. To gain a better understanding of the nar- 
rowing mechanism, many numerical and analytical stud- 
ies have been performed . From these studies it has 
emerged that if there are no correlations between the in- 
dividual environments of the different monomers, in the 
case of a one-dimensional arrangement of N monomers 
the line shape narrowing is usually of the order of 1/ y/~N 
compared to the width of the line shape of the non- 
interacting monomers and occurs if the coupling between 
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the monomers is much larger than the coupling to the en- 
vironment. However, recently it has emerg ed @,[TM1 
that for certain disorder distributions a totally different 
behaviour can be found. For example in the case of a 
Lorentzian distribution (i.e. a Lorentzian absorption line 
shape) there is no narrowing at all [f| [l^, [H| and for dis- 
tributions, which fall off slower than a Lorentzian, even 
a broadening can occur [l2j. 

In the present paper we will not consider static disor- 
der, but focus on the coupling of the excitation to vibra- 
tional modes of the environment. We use a commonly 
adopted model in which the excitation of a monomer 
is linear ly coup led to harmonic modes of its environ- 
ment 0, ll3Hl7l |. This model, restricted to only a single 
environmental mode, is often used to describe the in- 
fluence of a high energy vibrational mode in aggregates 
of organic dye molecules [H M, El- 

For one mode 

the resulting absorption spectrum of non-interacting 
monomers consists of delta-functions spaced by the fre- 
quency fl of the mode and the weight of the delta- 
functions, i.e. their "height", follow a Poisson distribu- 
tion 0, 0, [l(| whose standard deviation (width) is de- 
termined by the coupling strength between the electronic 
excitation and the vibrational mode (note that a Poisson 
distribution is completely determined by the standard 
deviation). For a chain of N monomers one then finds 
that in the limit of very strong interaction (i.e. the in- 
teraction between the monomers is much larger than the 
interaction between electronic excitation and vibrations) 
the absorption spectrum is again a Poissonian consist- 
ing of delta- functions spaced by ft, but with a standard 
deviation decreased by a factor 1/y/N. This leads to 
an increase of the intensity of the delta-function at low- 
est energy w.r.t. the delta-functions at higher energies. 
For more than one mode a similar pattern is found: the 
monomer line shape consists of a convolution of the Pois- 
sonians of the different modes, and consequently its vari- 
ance is given by the sum of the variances of the individual 
Poisson distributions. For the chain the standard devia- 
tion of each Poissonian is diminished by a factor 1 / y/N so 
that also the standard deviation of the whole peak struc- 
ture is diminished by a factor 1/y/N. When the density 
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of modes becomes large (continuous), the resulting ab- 
sorption spectrum becomes continuous, too. Often this 
spectrum is then dominated by a single band/peak which 
possesses a more or less pronounced shoulder. In such a 
case it is common to characterize the the width of this 
peak by its full width at half maximum (FWHM). In the 
following we will also use this measure. In the conclusion 
we will discuss the relation between different measures 
for the width of an absorption spectrum w.r.t. exchange 
narrowing. 

In this work we investigate the transition from the 
monomer line shape to the exchange narrowed line 
shape of the one-dimensional aggregate, considering a 
continuous spectral density of the environment, which 
leads to a continuous absorption lineshape of the 
monomer/ aggregate. We will focus in particular on the 
FWHM of the dominant absorption peak. We will find, 
that in analogy to the case of Lorentzian static disor- 
der also in the case of a Lorentzian monomer absorption 
line shape (corresponding to a Markovian environment) 
there is no exchange narrowing. On the contrary, in the 
case of a Lorentzian spectral density (which leads to an 
asymmetric structured monomer line shape) an enhanced 
narrowing is found. 

The paper is organized as follows: In Section [Til the 
Hamiltonian we will consider is introduced. In Sec- 
tion IIII1 the method to calculate the zero temperature 
absorption spectrum is briefly reviewed. The exchange 
narrowing is investigated in Section IIVI We conclude in 
Section [V] with a discussion of our findings. In the ap- 
pendix, an analytical formula for the monomer absorp- 
tion spectrum is given. 

II. THE MODEL HAMILTONIAN 

We consider an assembly of N two-level systems (i.e. 
one two- level system per monomer). The transition en- 
ergy between the ground state | </>^ ) and the excited state 
\4>n) °f monomer n is denoted by e n . We will consider 
absorption from the state where all monomers are in their 
ground state. The state in which monomer n is excited 
and all other monomers arc in their ground state is de- 
noted by 1 7r„ ) . In the following we will not consider 
states in which the whole system has more than one ex- 
citation. Expanding the Hamiltonian of the interacting 
monomers with respect to the "one-exciton" states 1 7r„ ) 
gives 

N 

+ V nm ) | TTn . ) ( 7T m |, (1) 

n,m— 1 

where the interaction V nm describes exchange of excita- 
tion between monomer n and m. Taking also the in- 
teraction with the environment into account the total 
Hamiltonian in the one-exciton subspace can be written 
as [l| 

H = H sys + Hi n t + H cnv . (2) 



The environment of each monomer is taken to be a set 
of harmonic modes described by the Hamiltonian 

N 

H cnv = y^y^ u n \a nX a n \. (3) 

n=l X 

Here a n \ denotes the annihilation operator of mode A of 
monomer n with frequency u> n \. Note that each monomer 
possesses its own environment which is independent of 
the environments of the other monomers. The coupling 
of an excitation on a monomer to the modes of its local 
environment is assumed to be linear and the correspond- 
ing interaction Hamiltonian is given by 

N 

Hint = | 7T„ ) ( 7T„ | ^ K nX (a' nX + 0,„ \ ) , (4) 

n=l A 

where the coupling constant k„\ specifies the strength 
of the coupling of the excitation on monomer n to the 
mode A with frequency ui n \ of the local environment. 
This interaction is conveniently described by introducing 
the spectral density [l4| 

J n {uj) = ^ \ K n\\ 2 S(uJ - LO n x) (5) 
A 

of monomer n. In this work we will consider J n (uj) to 
be a continuous function. The spectral density is closely 
related to the correlation function of the environment 
(bath correlation function). We restrict the discussion in 
this work to the case of zero temperature where the bath 
correlation function a„(r) reduces to 

a„(r) = J du J n (oj)e- iuT . (6) 

For a continuous spectral density J„(w) the bath corre- 
lation function is decaying to zero on a time scale de- 
termined by the width of J n (uj). When the environment 
has no memory, i.e. a„(r) cx 8(t), it is termed Marko- 
vian, otherwise it is non-Markovian. 

III. THE ABSORPTION SPECTRUM 

At zero temperature, initially the aggregate is in a 
state | g ) in which all the monomers are in their ground 
state | 0^ ) and also all environmental modes are in their 
ground state, i.e. 

N 

i.9>=i.9vib>ni^>> ( ? ) 

where | g v ih ) denotes the product of the vibrational 
ground states of all modes. For simplicity we take all 
transition dipoles of the monomers to be identical. Then 
the absorption strength for light with frequency v is given 
by Wm 

A{u) = Re / dt c wt ( * \e~ tHt I *o ) (8) 
Jo 
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with 



and 



*0) = IV>o)|<?vib) 



V>o) 



1 N 



(9) 



(10) 



The calculation of the aggregate absorption spectrum 
poses serious problems in the general case, due to the 
complex structure of the coupling to the environmental 
oscillators. In certain cases however, it is possible to 
solve the problem numerically or even analytically. In 
the following we will first consider the Markovian case, 
where wc show analytically that no narrowing occurs. 
To demonstrate the effect of anomalous strong narrow- 
ing we will concentrate on a spectral density which is a 
single Lorentzian. For such a spectral density we are able 
to calculate the absorption spectrum for quite large ag- 
gregates numerically exact, using the method described 
in Ref. S3. 



IV. INVESTIGATION OF THE EXCHANGE 
NARROWING 

To investigate the appearance of exchange narrow- 
ing we choose a one dimensional arrangement of the 
monomers and take the transition energies to be equal, 
e„ = e. Wc introduce periodic boundary conditions (ex- 
cept for the dimer) and take only coupling between near- 
est neighbors into account, which we assume to be equal 
for all neighboring monomers, i.e. T4i,n+i = V for all 
n. For simplicity (as already mentioned in the previous 
section) we also take all the transition dipole moments 
to be identical. For this arrangement, the absorption 
spectrum, in the limit of vanishing coupling to the envi- 
ronment, simply consists of one absorption line located 
at the band edge of the exciton-band pHE], [22j 



A{v) = 6(v-{e + C)), 



(11) 



where C = V in the case of a dimer and C — 2V for 
longer aggregates. 



A. Markovian case 

For a Markovian environment, one has a„(r) oc S(t), 
specifically we use a„(r) = 2T8(t). In this case, one can 



show [H Hi HH that ( * \e- lHt \ f ) needed in Eq. © 
is given by ( ip® \ tp(t) ) with | tp ) from Eq. (|TU|) and where 
the purely electronic state \ip(t)) can be obtained from 
solving 



with the initial condition \tp(t)} = \ipo). One can easily 
solve Eq. (fT2"j) and use Eq. ([5]) to obtain the absorption 
line shape 
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AM = 

^ {u-{e + C)) z +Y 2 



(13) 



This is a Lorentzian line shape centered at (e + V) in the 
case of a dimer and at (e + 2V) for N > 3. The width T is 
independent of the interaction V between the monomers 
and of the number N of the monomers, i.e. no narrowing 
occurs. A similar result has previously also been observed 
for the case of Lorentzian static disorder @, [1(1 EU ■ 



B. Lorentzian spectral density 

For the environment we consider a continuum of fre- 
quencies so that the spectral density J n (w) becomes a 
smooth function. We will consider a Lorentzian spectral 
density 



J n (oj) 



1 



-n 2 x- 



7 



(w - n) 2 + 7 2 



(14) 



d t \^(t)) = (-iH ays -T)\^(t)) 



(12) 



Note that in the limit 7 — ?► this spectral density 
becomes a delta function peaked at the frequency fi. 
This leads to the well investigated model, in which 
each monomer possesses one undamped vibrational mode 
of frequenc y f l that couples to the electronic excita- 
tion @, [lj, lH, GJ, [H Hi]. The quantity X appearing 
in Eq. (|14|) is the Huang-Rhys factor [26| that indicates 
the strength of the coupling between electronic excitation 
and the vibrational mode. The width 7 of the Lorentzian 
spectral density is the magnitude of the damping of the 
vibrational motion. 



1. The monomer spectrum 

For a Lorentzian spectral density the monomer absorp- 
tion line shape (i.e. for V = 0) can be calculated analyt- 
ically, as is done in the appendix. In Figure Q] the result- 
ing monomer spectrum is shown for various parameters 
X and 7. Here and throughout the paper the frequency 
f2, at which the Lorentzian spectral density is centered, 
is taken as the unit of energy. As the zero of energy 
we choose the monomer transition energy e, where the 
mean of the monomer spectrum is located. In the left 
column of Fig. [T]all spectra are for a width 7 = 0.1 and 
X increases from top to bottom. For increasing X the 
coupling of the electronic excitation to the environment 
becomes stronger resulting in the appearance of a vibra- 
tional progression. For a larger 7, i.e. a larger width of 
the spectral density, all peaks of the monomer spectrum 
become broader, as can be seen in the right column of 
Fig. [1] where 7 = 0.4. Note, that for small Huang-Rhys 
factors (X < 0.3) the monomer spectrum is dominated 
by one peak. 
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FIG. 1: Absorption spectrum of a monomer (i.e. V = 0) for various values of X (indicated on the figures) and for 7 = 0.1 (left 
column) and 7 = 0.4 (right column). The frequency Q, at which the Lorentzian spectral density is centered, is taken as the 
unit of energy (throughout the paper). 



2. The narrowing of the aggregate spectrum 

As an example, in the following investigations we will 
consider the case of the monomer spectrum Fig. QJ, ob- 
tained for X = 0.3 and 7 = 0.4. We now explore how this 
line shape narrows upon increasing the coupling strength 
\V\ between the monomers and examine its dependence 
on the number N of monomers. In Figure [5] aggregate 
spectra are shown for different N and for various values 
of the coupling V. For comparison, the corresponding 
monomer spectrum of Fig. [Tf is also shown (right most 
curve in each plot). For all spectra we chose V < (in 
molecular aggregates such a negative coupling leads to 
the formation of the famous J- band [22|). The larger the 
value of I V|, the more the spectrum is shifted to lower en- 
ergies. This is in accordance with a sum rule that states 
that the mean of the absorption line shape is shifted by 
C = 2V (and C = V for the dimer, i.e. N = 2, respec- 
tively) with respect to the mean of the monomer line 
shape [13, ■ From another sum rule it follows that 
the area of the absorption line shape is independent of 
V. Thus the increase of the height of the absorption 
peak with increasing \ V\ already indicates its narrowing. 
Going from the first to the third column of Fig. [2] the 
maximal \ V\ that is shown is increased, leading to a shift 
of the peak position to lower energies. 

To investigate the narrowing with increasing \ V\ in de- 
tail, consider first Fig. [2^,, where dimer spectra (N = 2) 
are shown for relatively small dipole-dipole interaction 
(1^1 < 2). One sees, that upon increasing \V\ the width 
of the absorption spectrum becomes gradually narrower. 
This narrowing takes place quite rapidly for \V\ < 1 and 



then slows down. For V ~ — 2 the spectrum has nearly 
obtained its asymptotic form and does not change ap- 
preciably upon increasing \V\ further. This can be seen 
clearly in Figs. and c, where dimer spectra for larger 
|V| are shown. 

In addition to the dimer case also spectra for N = 3 
(second row of Fig. [5]), N = 6 (third row) and N = 12 
(last row) are shown. These spectra are calculated for the 
same values of V as the corresponding dimer spectra of 
the first row. Since in a cyclic aggregate with N > 2 each 
monomer possesses two nearest neighbors, for the same V 
the effective coupling strength is now given by C = 2V 
leading to a shift of the mean of the absorption peak, 
that is twice as large as for the dimer [27l - l29j . Let us first 
consider the spectra for the case N — 3 (Figs.[5]i-f). Here 
a similar development of the width as for the dimer can be 
observed. The spectrum initially narrows very fast and 
reaches, for approximately V = —2, nearly its asymptotic 
line shape (as can be seen in Fig. [5^). However, the 
spectra are much narrower than the corresponding dimer 
spectra. 

When considering a longer aggregate with = 6 
(Figs. Hk-i) we find that one has to go to much larger 
I I than in the case N = 3 to reach the asymptotic line 
shape. For A^ = 6 it is reached roughly for V = — 12 (see 
Fig- HI). For even larger aggregates larger values of \V\ 
are needed to achieve the asymptotic line shape as one 
can see in Fig. for the case N — 12. Going from 
the top to the bottom row of Fig. [2] one observes that 
for the same V but increasing A^ the spectra also become 
significantly narrower. 

To investigate quantitatively the dependence of the 
narrowing of the spectrum on A^ and \V\, we numeri- 
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FIG. 2: Absorption spectra of different TV-mers (TV indicated on the figures) for various V; the corresponding monomer spectrum 
is that of Fig. QJ, which is shown here again for comparison (right most spectrum in each plot). In the first column the values 
of V are (from right to left): V = 0, —0.25, —0.5, —1, —2. In the second column in addition to the values of the first column 
also V — —3, —5 are shown. Third column: In addition spectra for V = —8, —12, —20 are shown. Note the different ranges of 
the abscissa of the individual plots. 



cally evaluated its full width at half maximum (FWHM). 
The results are shown in Figure [3J Here for different 
TV the FWHM is plotted against V. To show the nar- 
rowing for small |V| in closer detail, in the left column 
of Fig. [3] a smaller range of V is taken. While in the 
upper row the bare FWHM is plotted (in units of the 
FWHM of the monomer spectrum), the lower row shows 
the FWHM multiplied by the number of monomers, i.e. 
(TV-FWHM). In Fig. [3^, b one observes that the larger 
TV, the larger is the coupling strength \V\ needed for 
the FWHM curve to become horizontal, i.e. to reach the 
asymptotic FWHM. This is consistent with the observa- 
tion made by considering Figure O In Fig. one sees 
that the quantity (TV-FWHM) saturates for all TV at the 
FWHM of the monomer spectrum (which in Fig. [3J is 
taken as unity). That means that in the limit of strong 
interaction V, where the shape of the spectrum does not 
change anymore upon increasing \V\ further, the FWHM 
of the TV-mer spectrum is narrowed by the factor 1/TV 



w.r.t. the FWHM of the monomer spectrum. This is a 
much stronger narrowing than the narrowing by the fac- 
tor l/VTV found usually (the 1/a/TV narrowing has been 
found e.g. in the case of uncorrelated Gaussian disor- 
der for both the FWHM and the standard devia- 
tion of the absorption peak). However, it should be men- 
tioned that in contrast to this narrowing of the FWHM of 
the one dominant absorption peak by the factor 1/TV, as 
discussed above, the standard deviation of the spectrum 
only narrows by a factor 1/\N. A simple explanation 
for the narrowing of the FWHM will be given in the next 
section. 



C. Simple explanation of the anomalous strong 
narrowing 

From analytic arguments presented in Ref. @ one ex- 
pects that for | V| ^> A, where A denotes the width of the 
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FIG. 3: a) FWHM of JV-mer absorption spectrum (N in- 
dicated on the figures) over V; the corresponding monomer 
spectrum is that of Fig. [if and the FWHM of the N-mei 
spectrum is given in units of the FWHM of the monomer 
spectrum, c) (FWHM ■ N) over V. b, d) Same as a, c but for 
a larger range of V. The lines connecting the data points are 
only a guide for the eye. 

monomer absorption spectrum, the absorption spectrum 
of the interacting monomers (aggregate spectrum) has 
the same shape as the spectrum of the non-interacting 
monomers (monomer spectrum) obtained for a reduced 
coupling strength k\ = k\/^/N to the environmental 
modes. It follows from Eq. that according to the 
reduced coupling strength k\ one has an effective spec- 
tral density J{oS) = J(w)/N. This argumentation applies 
to arbitrary spectral densities, as long as the strong cou- 
pling condition \V\ S> A is fulfilled. 

For the Lorentzian spectral density Eq. ([T4]) that is 
considered here, this leads to a reduction of X by a fac- 
tor 1/N, i.e. X = X/N. Thus the aggregate spectrum 
in the limit of strong interaction V can be obtained by 
calculating the monomer spectrum for X = X/N and 
shifting it by C. That this reasoning is indeed correct is 
demonstrated in Figure [4] where for the case of strong 
interaction V, i.e. large |I^|/A, numerically calculated 
(exact) N-rner spectra (solid lines) are compared with 
the analytic monomer line shape obtained for X = X/N 
(dashed lines). The analytic monomer spectra are shifted 
by C = 2V (or C — V for the dimer respectively). In 
Figure H] three different values of V are chosen, namely 

V = — 5 (first column), V = — 20 (second column) and 

V = — 60 (third column). For the dimer the agree- 
ment between numerical and analytic spectrum is already 
nearly perfect for V = —20, as can be seen in Fig. HJd. 
However, for larger N one has to go to larger \V\ to 
achieve perfect agreement, shown for the cases N = 6 
and N = 12 in the third and fourth row of Fig.@] This is 
due to the fact observed already in the previous section 
that for larger N a larger coupling strength |V| is needed 



to reach the regime where the line shape has nearly con- 
verged, i.e. the strong coupling regime in which the ar- 
guments for the replacement X — > X/N of Ref. @ are 
valid. As described in Ref. the standard deviation of 
the monomer spectrum with X = X/N is reduced by 
a factor 1 / yN compared to the original monomer spec- 
trum with Huang- Rhys factor X . 

The narrowing of the FWHM of the aggregate spec- 
trum by a factor 1/N observed in Section [TV Bl can now 
be explained as follows: As is shown in the appendix, for 
the case X < 0.3 considered above, only one Lorentzian- 
like peak, whose width is proportional to X, dominates 
the monomer spectrum and contributes to the FWHM. 
Thus, when X is decreased by l/N, also the FWHM 
of the monomer spectrum decreases by 1/N. Since the 
shape of the aggregate spectrum in the limit of strong 
interaction V is that of the monomer spectrum with 
X —> X/N, the FWHM of the aggregate spectrum nar- 
rows by 1/N. 

V. CONCLUSIONS 

We have investigated the exchange narrowing for an 
one-dimensional aggregate interacting with Markovian 
and non-Markovian environments. In the case of a 
Markovian environment no narrowing occurred. For a 
non-Markovian environment with a Lorentzian spectral 
density we found enhanced narrowing with a factor 1 /N 
compared to the usual 1 / y/~N. This usual 1 / \J~N narrow- 
ing is found for uncorrelated Gaussian disorder (for the 
FWHM as well as for the standard deviation of the ab- 
sorption peak), but also for the standard deviation of the 
aggregate spectrum for the present model Hamiltonian of 
Section[Hj describing an exciton that couples to harmonic 
vibrations. Thus, some comments about the apparent 
contradiction between the 1/y/N narrowing, found e.g. 
in Ref. @, and the 1/N narrowing, found in the present 
work, are in order, since both rely on the same model 
Hamiltonian. It is important to note that the extent of 
the narrowing depends crucially on the definition of the 
"width" of the absorption spectrum. For example, from 
a sum rule it follows exactly that the standard deviation 
of the absorption spectrum of the aggregate is indepen- 
dent of the interaction V, in particular it is identical to 
the standard deviation of the monomer spectrum. Thus, 
no narrowing occurs, if one takes the standard deviation, 
which is a global quantity of the whole spectrum, as a 
measure of the width. It must be pointed out, that also 
weak absorption lines, well separated from the absorp- 
tion peak one is interested in, contribute to the standard 
deviation. In Ref. Q such small contributions very far 
away from the absorption peak considered are neglected. 
One then finds that the coupling strength to each vibra- 
tional mode is decreased by a factor 1/^/N, leading to 
an absorption line shape that is a Poissonian, but with 
a standard deviation decreased by 1/VN w.r.t. the Pois- 
sonian of the monomer spectrum. Thus, one obtains an 
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FIG. 4: N-mer spectra for large | V|/A, where A is the width of the monomer spectrum. Comparison of numerically calculated 
(exact) iV-mer spectra (solid lines) with the analytical monomer line shape obtained from the replacement X — > X/N and 
shifted by C (dashed). Left column: V = —5, middle column: V = —20, right column: V = —60. Note the different energy 
and absorption scales. 



overall 1 / y/N narrowing - again based on the standard 
deviation which gives large weight to small contributions 
far away from the mean. 

In the present paper we have considered the full width 
at half maximum (FWHM) of an absorption peak, which 
is a local measure of the width. It is applicable when- 
ever there is a well defined peak, as in the present study. 
By choosing this measure we neglect the wing at high 
energies, which can be seen in all spectra and which con- 
tributes the additional width when the standard devia- 
tion is considered. 

Furthermore, it should be mentioned that for larger 
Huang-Rhys factors, when more than one peak of the 
vibrational progression in the monomer spectrum con- 
tribute to the FWHM, the narrowing of the FWHM will 
not have the simple 1/N dependence anymore (see ap- 
pendix). 



While in this work we have focused on a Lorentzian 
spectral density, the main results, namely an enhanced 
narrowing compared to the 1 / case, are expected to 
remain valid also for other forms of the spectral density 
(as discussed in Section TlV Cp . We have performed cal- 
culations also for an ohmic spectral density with expo- 
nential cutoff and found similar results as discussed here 
for the Lorentzian spectral density. 

The findings of this work demonstrate, that one should 
be careful when deducing quantities such as the number 
of coherently interacting monomers from the narrowing 
of the absorption spectrum, since the exact value of the 
narrowing strongly depends not only on the model as- 
sumed (e.g. static disorder, single vibrations, a continu- 
ous spectral density, etc.) but also on the definition of 
the width. 
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Appendix A: Analytical formula for the monomer 
absorption spectrum 

The monomer absorption spectrum is given by (cf. 
Eq. ©) 

A{u) = Re / dt e wt c(t), (Al) 
Jo 

where the correlation function c(t) = ( "Jo \c~ lHt | "Jo ) can 
be determined analytically for various concrete bath cor- 
relation functions a(r) = J dtu e~ lUT J(u>) (here for zero 
temperature; J{uj) is the spectral density of the environ- 
mental modes). In the case of the Lorcntzian spectral 
density Eq. (|14[) . that we consider in this work, we have 

a(r) = T e - mT -~< T (A2) 

with r = Q 2 X (note that since in this work we always 
take = 1, the parameter T is equal to the Huang- 
Rhys-factor X). For this a(r) the correlation function 
c{t) needed in Eq. (|A1|) can be evaluated analytically to 
give [H 




where e is the transition energy of the monomer and with 
w = iO + 7. Inserting Eq. (|A3[) into Eq. (|A1[) leads to an 
analytic expression for the absorption spectrum 

AM = e^yr g MO ( _ 1 sign(7 2_ fl2) ) \ 

k=0 ' ' 

(A4) 

with the definitions 

7 = r 7 /( 7 2 + 2 ) (A5) 

n = ra/( 7 2 + o 2 ) (A6) 

and with 

Ak{v) = cos(# — kq) Im L^(v) + sm(9 — kq) Re L/.(u), 

(A7) 



where 9 = 2 7 fi/r and q = arctan(— 27O/ (7 2 — f2 2 )). Here 

L k M = — ^ (A8) 

is a complex Lorentzian line shape function. The spec- 
trum Eq. (|A4I) consists of a sum of terms where the line 
shape of each term is determined by Lf.{y), The imagi- 
nary part of L^(y) is a Lorcntzian centered at (fcO— fl+e) 
and with a width (7 + k"f). Thus the Lorentzians for the 
different k arc spaced with a distance O and their width 
increases with increasing k. The real part of L^(y) is 
given by the same Lorentzian but multiplied with the 
linear function (v — k£l + (l — e)/(7 + kj). 

Since the first peak (k = 0) is of particular interest in 
this paper (because it dominates the spectrum for small 
T) we discuss it in more detail. Denoting the correspond- 
ing term in the sum in Eq. (|A4[) by Aq(v) we have 

AoM = f 7 ^^-2 fcos(g)+sin(g) ' / + "~ g > ) . 

(u + O - e) 2 + 7 2 y 7 J 

(A9) 

For the values taken in the investigations of Section HVl 
i.e. r < 0.3 and 7 = 0.4, the line shape of this peak is to a 
good approximation a Lorentzian with a width (HWHM) 
given by 7. Since 7 oc T, replacing T by T/N to obtain 
the line shape of the iV-mer in the strong coupling limit 
(as has been done in Section HV C[) one finds a narrowing 
of the 7V-mer peak by a factor 1/N w.r.t. the monomer. 
However, for a stronger coupling T to the environment, 
when beside the k = peak also other peaks contribute to 
the FWHM of the spectrum, the narrowing of the FWHM 
does not have the simple 1/N dependence anymore. 
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